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Abstract 

This paper is concerned with front-like entire solutions for monostable reaction- 

■ diffusion systems with cooperative and non-cooperative nonlinearities. In the cooper- 
ed . ative case, the existence and asymptotic behavior of spatially independent solutions 

(SIS) are hrst proved. Combining a SIS and traveling fronts with different wave speeds 

■ and directions, the existence and various qualitative properties of entire solutions are 
then established using comparison principle. In the non-cooperative case, we intro- 
duce two auxiliary cooperative systems and establish some comparison arguments for 

■ the three systems. The existence of entire solutions is then proved via the traveling 
fronts and SIS of the auxiliary systems. Our results are applied to some biological and 
epidemiological models. To the best of our knowledge, it is the first work to study the 

■ entire solutions of non-cooperative reaction-diffusion systems. 
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1 Introduction 

This paper is concerned with entire solutions of the following m-dimensional reaction- 
diffusion system in M^: 

u t = DAu + f(u), xeR N ,teR, (1.1) 

where m, iV € N, 

u= {ui,--- ,u m ), / = (/!,••• ,/ m ), D = diag(di, • • • ,d m ), 

and (di, ■ ■ ■ , d m ) S> := (0, • • • , 0) € M m . Here and in what follows, we always use the 
usual notations for the standard ordering in W 11 . As usual, system (jl.ip is said to be 
cooperative on / C W 71 if each fi(u) is non-decreasing in Uj on / for 1 < j ^ i < m; 
otherwise, it is said to be non-cooperative on /. 

One important topic for reaction-diffusion systems is the traveling wave solution that 
describes the phenomenon of wave propagation. In the past decades, many studies have 
led to almost complete description of traveling wave solutions of (jl.ip with cooperative 
nonlinearity [2H Ell EH SHIS]. For example, Volpert et al. [M] gave a complete result 
about the monostable and bistable traveling fronts, and Tsai [31J investigated the global 
exponential stability of the bistable traveling fronts. In a series of papers, Weinberger, 
Lewis and Li (2TJUTH32] studied spreading speeds and traveling fronts for general cooper- 
ative recursion systems. Related results on scalar non-monotone evolution equations, we 
refer to [HIZllSSlllSlESlESlllHrllS] . 

In addition to traveling wave solutions, another important topic in diffusion systems is 
the interactions of them, which is crucially related to the pattern formation problem. We 
refer to [5,6,18,26] for more details. Mathematically, this phenomenon can be described 
by the so-called front-like entire solution that is defined for all space and time and behaves 
like a combination of traveling fronts as t — > — oo. On the other hand, from the dynamical 
points of view, the study of entire solutions is essential for a full understanding of the 
transient dynamics and the structures of the global attractor \27\ . In the recent years, 
there were many works devoted to the interactions of traveling fronts and entire solutions 
for scalar reaction-diffusion (both spatially continuous and discrete) equations with and 
without delays, see e.g., [ZHillllllllll^^ 

More recently, Morita and Tachibana [28j, Guo and Wu [H], and Wang and Lv [37] and 
Wu [36] extended the existence of entire solutions for scalar equations to some specific two 
component cooperative reaction-diffusion model systems. The basic idea in these studies, 
similar to [SJHSIEHIEH] , is to use traveling fronts propagating from both directions of the 
x-axis to build sub- and supersolutions, and then prove the existence results by employing 
comparison principle. Unfortunately, it seems difficult, if not impossible, to construct such 
supersolutions for the m- component reaction-diffusion system (jl.lj) . In fact, to the best 
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of our knowledge, there has been no results on the entire solutions for general cooperative 
reaction-diffusion systems and non-cooperative systems. 

The purpose of the paper is to consider the entire solutions of system (jl.ip with 
cooperative or non-cooperative nonlinearity. In the cooperative case, the existence and 
asymptotic behavior of spatially independent solutions are first proved. Since it is difficult 
to use traveling fronts to construct supersolutions for the general m-component system, 
we extend the method developed in [15] for scalar KPP equations to system (jl.ip . More 
precisely, we construct appropriate upper estimates by virtue of the exact asymptotic 
behavior of the traveling fronts and spatially independent solution, and then prove the 
existence and qualitative features of entire solutions using comparison principle (Theorems 
l2.9l and l2.lUJ) . Although the method is inspired by the work of Hamel and Nadirashvili [15J, 
the technical details are different. In [15] . the upper estimates were proved by the solution 
formulation of the linearization of the scalar KPP equation at the trivial equilibrium. 
Contrasting to [15], we use a general comparison principle to prove the upper estimates 
(Lemma I2.12j) . Recently, the method was successfully applied in our previous work [38] 
to a multi-type SIS nonlocal epidemic model. 

For the non-cooperative reaction-diffusion systems, we introduce two auxiliary coop- 
erative systems, one lies above and another below of system (jl.ip . which were used by 
Wang [36] and several references therein to obtain the existence of traveling wave solu- 
tions, and establish some comparison arguments for the three systems. Combining the 
traveling fronts and spatially independent solution of the lower system and their exact 
asymptotic behavior, we then build appropriate subsolutions and upper estimates of the 
auxiliary lower and upper systems using the comparison theorem, respectively, and prove 
the existence and qualitative properties of entire solutions of (II. ip with non-cooperative 
nonlinearity (Theorem 13. 6j) , To the best of our knowledge, it is the first work to study the 
entire solutions of non-cooperative reaction-diffusion systems. 

In biology and epidemiology, there are quite a few reaction-diffusion model systems of 
the form (jl.ip with cooperative or non-cooperative nonlinearities. We shall illustrate our 
main results by discussing the following models in [T, 20, 32, 331136]. 

A. A Buffered System. In |32jl33j. Tsai and Sneyd presented a buffered system: 

f d t ui = dAui + g( Ul ) + IX#rO° - v i) ~ k t u i v i]i ^ 2 ) 

| d t Vi = diAvi + k^(b® - vi) - kfu\Vi, i = 1, • ■ • , n, 

where d, k^ , > and d{ > are given parameters. They studied the existence, unique- 
ness and stability of traveling fronts of (jl.2p by taking the typical bistable nonlinearity 
for the function g, i.e. g{u\) = ui(ui — a)(l — ui) for some a € (0, 1). Note that (jl.2h can 
be transformed to a cooperative system on R + x Hr=i[0> under the change of variable 
Ui = Vi, i = 1, • • • , n. Other results related to the buffered system, we refer to [7jll3pi9] 
and the references therein. 
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B. An Epidemic Model. To study the fecally-orally transmitted diseases in the 
European Mediterranean regions, Capasso and Maddalena [TJ introduced the epidemic 
model: 

d t ui = diAui - auiii + a 12 u 2 , ^ ^ 

d t u 2 = d 2 Au 2 - a 22 u 2 + g{ui), 

where d%, a\x, a\ 2 , a>22 > and d 2 > are given parameters. The function g{u\) decribes 
the infection rate of human under the assumption that total susceptible human population 
is constant. In general, g(-) is increasing on [0,+oo). But, if the "psychological" effect is 
considered (see, e.g., Xiao and Ruan [39]), then g(-) is a unimodal curve on [0, +oo), that is, 
g(-) achieves its maximum at some u max > 0, and is increasing on [0, it max ] and decreasing 
on [« max , +oo). When d 2 = and g is monotone, Xu and Zhao |50j proved the existence, 
uniqueness and stability of bistable traveling fronts of (jl.3p and Zhao and Wang [53J 
established the existence and non-existence of monostable traveling fronts. These results 
were then extended by Wu and Liu [44] to the non-monotone case by constructing two 
auxiliary monotone integral equations. 

C. A Population Model. Weinberger, Kawasaki and Shigesada [20] discussed the 
reaction-diffusion model which describes the interaction between ungulates with linear 
density u\ and grass with linear density u 2 : 

d t ui = diAui + ui[-a - 5ui + riu 2 ], 
d t u 2 = d 2 Au 2 + r 2 u 2 [l - u 2 + h(ui)}, 

where d\, d 2 , r%, r 2 , a, 5 are all positive parameters. The function h(u\) models the in- 
crease in the specific growth rate of the grass due to the presence of ungulates. When the 
density u\ is small the net effect of ungulates is increasingly beneficial, but as the density 
increases above a certain value, the benefits decrease with increasing. In [20], Weinberger, 
Kawasaki and Shigesada established the spreading speeds for f|l .4[) by employing com- 
parison methods. Taking the non-monotone Ricker function uie~~ Ul as h(u\), Wang [36] 
further characterized the spreading speed as the slowest speed of traveling wave solutions. 
Throughout this paper, we always make the following assumptions: 

(A ) There exists K » such that /(0) = /(K) = 0, / € C 2 ([0, K], R m ) and there is no 
other positive equilibrium of / between and K. 

(Ai) One of the following holds: 

(a) The matrix /'(0) is cooperative and irreducible with s(/'(0)) > 0, where 

s(/'(0)) := max{KA : det(AJ - /'(0)) = 0}; 

(b) For each A > 0, A(X) := DX 2 + /'(0) is in block lower triangular form, the 

first diagonal block has a positive principal eigenvalue M(A), and M(A) is 
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strictly larger than the principal eigenvalues of all other diagonal blocks. In 
addition, there is a positive eigenvector v(X) = (v\(X), ■ ■ ■ ,v m (X)) ^> of A(X) 
corresponding to M(X) and v(X) is continuous with respect to A. 

We mention that a square matrix is called to be cooperative if all off-diagonal entries 
are non-negative, and irreducible if it cannot be placed into block lower-triangular form 
by simultaneous row/column permutations (Smith |29j). 

If (Ai)(b) holds, by the argument of |36j Lemma 1.1], there exist two numbers c* > 
and A* > such that 

c , _ mim _ inf m, (1 . 5) 

A* a>o A 

and for any c > c*, there exists Ai := Ai(c) G (0, A*) such that M(Ai) = cAi and 
M(A) < cA for any A G (Ai,A*]. 

If (Ai)(a) holds, then the matrix A(X) = -DA 2 +/'(0) is also cooperative and irreducible. 
Hence 

M(A) = s{A(X)) := max{KA : det(AJ - A(X)) = 0} 

is a simple eigenvalue of A(X) with an eigenvector v(X) = (vi(X),--- ,v m (X)) » 0. In 
addition, M(A) = s(A(X)) > a(/'(0)) > for any A > (see e.g., [221 Corollary 4.3.2]). 
From the argument of Lemma 2.1], there also exist c* > and A* > such that (11.51) 
holds, and for any c > c*, there exists Ai := Ai(c) G (0, A*) such that M(Ai) = cAi and 
M(A) < cA for any A G (Ai,A*]. 

The rest of the paper is organized as follows. In Section 2, we consider the entire 
solutions of system (jl.lj) with monostable and cooperative nonlinearity (Theorems 12.91 
and I2.10p . Section 3 is devoted to the entire solutions of (II. ip with monostable and non- 
cooperative nonlinearity (Theorem l3.6p . In Section 4, we apply our abstract results to the 
above models (|1.2jl - (|1.4j) . Finally, conclusions and discussions are given in Section 5. 



2 Entire solutions for cooperative systems 

In this section, we consider the entire solutions of (jl.ip with monostable and cooper- 
ative nonlinearity. In addition to (Ao) and (Ai), we also need the following assumptions: 

(A2) System (jl.lj) is cooperative on [0,K], that is, djfi(u) > for all u G [0,K] and 
1 < j ^ i < m. 

(A 3 ) For any k G Z + , pi,-- - , pu > and Ai,-- - , X k G [0,A*], 

/(min{K, pn;(Ai) + • • • + p k v(X k )}) < /'(0) [piv(Xi) + ■■■ + p k v(X k )) . 



Here, v(X) ^> is the eigenvector of A(X) corresponding to M (A). 
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Remark 2.1 It is easily seen that if f(u) < f'(0)u for u G [0, K], then (A3) holds 
spontaneously. We also note that if / is defined on [0, +oo) m , then (A3) can be replaced 
by (A 3 )*: 

(A 3 )* For any k G Z+, pi, ■ ■ ■ , p k > and A x , • • • ,A fe G [0,A*], 

/(pi«(Ai) + • • • + Pfcv(Afe)) < /'(0) [pi«(Ai) + • • • + Pkv(Xk)] ■ 

From the arguments of Theorem 3.1] and [36, Theorem 2.1], we have the following 
result. 

Proposition 2.2 Let (Ao)-(As) hold. For every c > c* and v G M N with \\is\\ = 1, (II. ID 
admits a traveling front 

*c(0 = (</>l,c(0> • • • > 0m,c(O)> ^ = X • 1/ + Ci, 

which satisfies & c (— 00) = 0, <3? c (+oo) = K and «$<;(•) ^ 0. Furthermore, there holds 

lim $ c (£)e~ Al(c)§ = u(Ai(c)) and $ c (f) < t>(Ai(c))e Al(c)f /or aZZ ^1. 

£^—00 

In the remainder of this section, we first give some comparison theorems for sub and 
supersolutions of (|l.lj) . We then state the main results for the cooperative system (Theo- 
rems 12.91 and I2.10D and establish the existence and asymptotic behavior of spatially inde- 
pendent solutions. Finally, we prove Theorems 12.91 and 12.101 by constructing appropriate 
subsolutions and upper estimates and using a general comparison principle. 

2.1 Preliminaries 

Consider the initial value problem of (jl.ip with initial condition: 

u(x,t) = ip(x), x<ER N , (2.1) 

where r G M is an any given constant. 

Let X = BUC(R" ZV , M m ) be the Banach space of all bounded and uniformly con- 
tinuous functions from M. N into W 71 with the supremum norm || • \\x- For simplicity, 
we denote W = [0,K] and [0,K] X = {(f> G X : < <j)(x) < K, x G R N }. Take L = 
{\difi(u)\\u G [0,K]} and define 

Q{u) = (Q x (u), --- , Q m {u)) = f(u) + Lu, u G W. 

Clearly, Q(u) is non-decreasing in u for u (zW. We further define a family of linear operator 



T(t) = diag(Ti(t), • • • ,T m (t)) :X^X,t>0, 



(2.2) 
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by Ti(0) = I and 

(Ti(t)<f>)(x) = e~ Lt [ ^i{y,t)(f>(x - y)dy, Vx t> 0, <j>(x) € BUCQR 

Jr n 



where 



1 



■ cxp 



1, • • • ,m. 



(4di7rt) N / 2 I 4di* 
The definitions of sub- and supersolutions of (jl.ip are given as follows. 

Definition 2.3 A continuous function u = (u±, • • • , u m ) : ~R N x [t, +oo) — > W is called a 
super 'solution of ( fi. 1}) on [r, +oo) if 

u(x,t)>T(t-T)u{x,r)+ [ T(t- s)Q{u(x,s))ds, VxeR N ,t>r, (2.3) 



A subsolution of M.l}) is defined by reversing the inequality. 

Remark 2.4 Let w = (w\, • • • , w m ) : R N x [r, +oo) — > W be a continuous function with 



the property that wi is C 1 in t and C 2 in x. It is easy to see that if w satisfies 

w t > {or <)DAw + f{w), Vx eR N ,t> t, 
then w is a supersolution (or subsolution) of on [r, +oo). 

By Definition 12. 3| we have the following results, see e.g., Fang and Zhao [7J. 

Lemma 2.5 (i) For any cp £ [0,K]x, (COP admits an unique classical solution u(x,t;(p) 
satisfying u(x, r; ip) = ip(x) and < u(x, t;ip) <~K for all x 6 1^ and t > r . 
(ii) Let w + (x,t) and w~(x,t) be a supersolution and a subsolution of respectively. 
If w + (-,t) > u>~(-,r), then w + (-,t) > w~(-,t) for all t > r. 

The following result follows from the standard parabolic estimates (Friedman |10|). see 
also Wang et al. [38, Proposition 4.3]. 

Lemma 2.6 Suppose that u(x, t; ip) is a solution of with the initial value ip € 

[0, K]x, then there exists a positive constant Mi, independent of r and (p, such that for 
any x € R N andt>r + l, 



du . . 


< Mi, 


d 2 u , , 
dtxfi^ 


< Mi, 


du . 


<Mi, 


d 2 u . , 
dx/^ 


<Mi, 



1i 



dt 2{ 

d 2 u 



dx i X j 



x,t;<p) 
-(x,t;tp) 



<Mi, 
< Mi 



<9 3 -u , , 



<Mi, 



d 3 i 



dx 2 xi 



{x,t;<p) 



<Mi, Vi,j = !,-•• ,iV. 
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Similar to Lemma l2.5f ii). we have the following result. 

Lemma 2.7 Let u + G C(R N x [t,+oo), [0, +oo) m ) and 

u~ G C(R N x [t, +oo), (-oo, Ki] x • • • x (-oo, K m }) 

be such that u + (-,t) > u~(-,t) and 

uf > DAu + + /' (0)u + , Vx G R N , t > r, 
Ut <DAu~ + f'{0)u~, Vx<ER N ,t>r. 

Then, u + (x,t) > u~(x,t) for all x G R N and t > r . 

2.2 Main results for cooperative systems 

Before to state our main results, we give the following definition and notation. 

Definition 2.8 Let n G N and p,po G R n . We say that the functions W p (x,t) = 
(Wi ;p (x,t), ■ ■ ■ ,W m - p (x,t)) converge to W Po (x,t) = (Wi ;Po (x, t), ■ ■ ■ ,W m - Po (x,t)) as p -4 
Pq in the sense of the topology T if, for any compact set S C R^"^ 1 , the functions W p , 
dtW p , d Xi W p , d x 2W p , i = I,--- ,N, converge uniformly in S to W po , d t W P0 , d Xi W po , 
d x 2W P0 , i = !,-■■ , N, as p ^ p . 

Notation: For any I G Z + , Vi G R , i = 1, • • • ,1, A G R and a G R, denote the regions 
T\a andTi a ,i = l,... ,1 + 1, by 

T A a := {x G R N \x ■ Vi > A] x [a, +oo), % = 1, ■ ■ ■ ,1, T^ 1 := R N x [a, +oo), 
f\ a := {x G R N \x -Ui<A}x (-oo,a], i = !,-■■ ,1, f 1 ^ := R N x (-00,0]. 

Now, we state the main results for the cooperative system as follows. 

Theorem 2.9 Let (Ao)-(A3) hold. Then, for any I G Z+, v x ,--- ,vi G R N with \\ui\\ = I, 

hi,-- - ,hi +1 G R, ci, • • • , q > c* ; and Xl, • • ■ ,Xm G {0, 1} wt/i XH l"Xz+i > 2, i/iere 

exists an entire solution U p (x,t) := (Ui ;p (x,t), ■ ■ ■ ,U m - p (x,t)) of (jl.ip such that 

u(x,t) <U p (x,t) < mm {K,U(x,t)}, V(i,i)eR w , (2.4) 

where p := p Xl ,-, Xl+1 = (xici, Xihi, Xi"i, • • • ,XicuXih,Xi^uXi+ih+i) and 

u(x,t) := max { . max Yi$cj (z ■ + qi + hi),xi+iT(t + /ty+i)}, 
1 

U(x,t) := J2xMXi(ci))e Xl ^ x ^ +c ' t+h ^ + X i+iv* e x *^ t+h '+^ . 
i=l 

Here, T(t) is the spatially independent solution of (11. ip decided in Lemma \2.11\ A* = M(0) 
and v* = v(0). 

Furthermore, the following statements hold: 
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(i) « U p (x,t) < K and ^U p (x,t) » for all (x,t) € R N+1 . 

(ii) limt_ > _ 0O sup|| a ,||< vl = for any AgN. 

(iii) If Xl+1 = °> then lim t ^ +00 sup^y^ \\U p (x, t) - K|| =0 for any A £ R + , and if 

Xi+i = 1, then limj^+oo sup^gRiv ||f7 p (a;,i) -K|| = 0. 

(iv) If xi+i = 1; then for every x £ M> N , 

U p (x, t) ~ r(t + ftj+i) ~ ^ e A *(*+^+i) as t _> _oo. 

(v) If xi+i = 0; then for every x £ M^, 

as i —> — oo, 

w/iere 0(ci, • • • ,q) = min (ciAi(ci), • • • , qAi(q)}. 

(vi) For any (x,t) € M Ar+1 , U p (x,t) is increasing with respect to hi, i = 1, • • • , I + 1. 

(vii) When N = 2 and I = 2, let us denote vi = (cos 0i, sin QA, 6i G [0, 2-7r), i = 1,2. Tjf 
Xi = 0, i/ien 

9 rr ^_ / >0, 2 € [0,f]U[f ,27r]; 



and 



JLu( xt )-l ^ e M; 

dx 2 pl ' j ~\ <0, 2 G[7r,27r]. 

Similar results hold true for N = I = 2 and \2 = 0. 

(viii) U p (x,t) converges to K as /ij — >• +oo m T and uniformly on (x,t) £ 1^ /or any 
A,a£l,i = V" ,1 + 1. 

According to the assumption xi, • • • , Xz+i £ {0j 1} with xi + ' ' ' + Xi+i > 2 in Theorem 
12.91 we denote the entire solution U p (x, t) of (jl.ip by 

U po (x,t), if (xi,-- - ,Xi+l) = (V" j 1 ); 
U p (x,t) := < U p .(x,t), if (xir-- ,Xz+i) = (!,••• ,1,0*, l,--- , l), i = 1, • • • ,Z; 
f^Pi + i(M) 3 if (xir-- ,Xi+i) = C 1 , • ■ • >l,0j+i), 

(2.5) 

where po = (ci, hi, z'i, • • • ,c h h h u h hi^, pi +1 = (ci, /ti, z/i, • • • ,q,h z ,i^,0) and 

K = (ci,/ii,z^i, • • • ,Ci__i,/ii_i,fi_i,0,0,0,c i+ i,/ii + i,^ + i,--- ,ci,hi,u h hi + i), % = !,■■■ ,1. 
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Moreover, we denote 



U p (x,t) := < 



where 



U Pij (x,t), if (xi,-- - ,Xl+i) = (1>" - , 1, 0i, 1, • • • ,l,0j,l,-- - ,1,1), 

i<»,i<i; 

tfpi.i+iOM), if (xi,--- ,Xl+i) = (!,••■ , 1, i; 1, • • • ,l,0;+i), i = l,--- ,1, 

(2.6) 



— (ci, hi, v\, ■ ■ ■ , Cj_i, /ij-i, fj-i, 0, 0, 0, Cj+i, /ij+i, i^i+i, • • • , 

c j-i,hj-i,Vj-i,0,0,0,Cj +1 ,hj + i,v j+ i, - ■ ■ ,ci,hi,ui,hi +1 ), 1 < i < j < I, 

and 

Pi,l+l = {ci,hx,vi, ■ ■ ■ ,Ci_i,/ij_i,fj_i,0, 0,0,c i+ i,h i+ i,v i+ i, ■ ■ ■ ,ci,hi,vi,0), i = 1, ■ ■ ■ ,1. 

Then we have the following convergence results. 

Theorem 2.10 Assume (Ao)-(A3). Assume further that f'(u) < /'(0) for u € [0,K]. 
Then, from (j2.5[) and (|2.6p . the following properties hold. 

(i) For any A G Z and a € R, U po (x,t) converges to 

U Pi (x, t) as hi — > — oo in T, and uniformly on (x, i) £ , i = 1, ■ • ■ , i + 1. 

(ii) For any AgR and a € M, U Pi (x,t)(i = 1, • • • , /) converges to 

U Pij (x,t) as hj — > — oo in T, and uniformly on (x,t) € T^ a , 1 < i 7^ j < I.; 
U Pil+1 {x,t) as hi + i — > —00 in T, and uniformly on (x,t) € T 1 ^. 

(iii) For any A £ M and a£l, U Pl+1 (x,t) converges to 

U Pi i+1 (x, t) as /ij — >■ —00 in T, and uniformly on (x, t) E T\ a , i = 1, • • ■ , /. 

(iv) For any hi, - ■ ■ ,h\,h\,--- ,h^ £ R, t/iere exists {xq, to) £ M Ar+1 , depending on c\, - ■ ■ , 

q, h\, - ■ ■ ,hi,h\,--- ,h* t such that 

U Pt+1 (x,t) = U pt+i (x + x ,t + t ) for all (x, t) £ R N+1 . 
Here, p* l+1 = (c x , h\, i>i, ■ ■ ■ , ci,h*,vi, 0) . 



Entire solutions in monostable R-D systems 



11 



2.3 Existence of spatially independent solutions 

In this subsection, we consider the spatially independent solutions of (|l.ip connecting 
and K, that is, solutions of the following ordinary differential problem: 

^ = /(r(t)),tei, (2.7) 

T(-oo) = 0, T(+oo) = K, (2.8) 

where r = (Ti, • • • ,T m ) and / = (ft, ■ ■ ■ , f m ). Recall that W = [0, K]. 

Note that (|2.7p is a cooperative and irreducible system. The existence of such a 
heteroclinic orbit T(t) can be established by using the theory of monotone dynamical 
systems (see Smith [29] and Zhao [52]). However, these results do not give the exponential 
decay rate of the solution at minus infinity. To overcome the shortcoming, we shall use the 
standard technique of monotone iteration scheme to prove the existence and asymptotic 
behavior of the solutions of (12.71) and (12. 81). 



Lemma 2.11 Let (A )-(A 3 ) hold. There exists a solution T(t) : R -> W of ([277)) and 
fl%lD such that 

r'(i) > 0, lim r(t)e~ A ** = v* , and T(t) < e A *V for all t G R, 
where X* = M(0) and v* = v(0). 

Proof. Since the method is standard, we only sketch the outline. Let C(R,R m ) be the 
spaces of continuous vector- valued functions on R. Define the operator F = (F\, • • • , F m ) : 
C(R,W) -> C(R,R m ) by 

Fi(u)(t)= [ e- L{t - s) Qi{u){s)ds, i = 1,--- ,m. 

J — oo 

Recall that 

L= max {\difi(u)\\u G W} and Qi(u)(t) = fi(u(t)) + Liii(t), i = l,---,m. 

It is easy to verify that each Qi(-) is a nondecreasing map form C(R, W) to C(R,R) with 
respect to the point-wise ordering. The remainder of the proof is divided into the following 
there steps. 

Step 1. The following observation is straightforward. 

(i) F : C(R,W) -> C(R,W); 

(ii) F((j))(t) > F(ip)(t) for </>,ip G C(R,W) with > ^(i); 

(iii) F(cj))(t) is increasing in R for G C(R, W) with is increasing in R. 
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Step 2. For any fixed e G (l, 2) and sufficiently large q > 1, define two functions as follows: 
^(t) = (&(*),■■■ A m {t)) and^(t) = (*),••• ,£jt)), 

where 

= min {Xi, <e A **} and 0.(i) = max jo, <e A ** - qv*e eXH } ,t€R. 

Then, by direct computations, we obtain 

< < 0(t) < K, F($)(t) < 4>{t) and > (jAt) for all i € R. 

S'tep 5. Using the monotone iteration technique, we can show that equation (|2,7|) admits 
a solution T(t) which satisfies 

T'(t) > and 0(i) < r(t) < 0(t) for all t G K. 

Thus, r(-oo) = 0, T(+oo) G (0,K] and 

lim r(t)e' A ** = v*, < r(t) < e A *V for all f G R. 

t— >— oo 

Moreover, one can easily verify that r(+oo) = K for all t G R. 

Next, we show that T'(t) > for all t G R. Since djfi(u) > for all u G [0,K] and 
1 < j 7^ i < ™, by (|2.7p . we have 

r?(t) = d 1 f i (T(t))T' 1 (t) + --- + d m f i (r(t))r' m (t) 

> m r-(t), Vt G R, 
where mo = min u G W\. Thus, for any r G R, we obtain 

i=l, ■•■ ,m 

r^(t)>r^(T)e m °^, Vt>r, i = l,-.. ,m. (2.9) 

Suppose for the contrary that there exist io G {1, ■ ■ ■ , m} and to G R such that r^ o (io) = 0, 
it then follows from that r- (r) = for all r < i . Thus, r io (r) = r io (i ) for all 
t < to and hence < rj (io) = Fj (— oo) = 0. This contradiction shows that T'(t) 3> 
for all t€K. The proof is complete. □ 

2.4 Proofs of Theorems and I2TTU1 

In this subsection, we will use the results of previous subsections to obtain an appro- 
priate upper estimate for solutions of (|1.1|) and then prove Theorems 12.91 and I2.1U1 
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For any l,n £ Z + , v\, ■ ■ ■ ,Vi € with ||z^|| = 1, hi, ■ ■ ■ , hi + i £ 1, q, • ■ • , q > c,, 
and xi, ■ ■ ■ , Xi+i € {0, 1} with xi H h Xz+i > 2, we denote 

(p n (x) := max 



/ max^^c^ • Vi - Cin + /tj) , xi+iT(-n + fy+i) j , 
u(x,t) :=maxi max Xi^ Ci (x ■ Vi + <%t + hA , Xl+iT{t + \ , t> 



-n. 



Let U n (x, t) = (Ui(x, t), ■ ■ ■ , Um(x, t)) be the unique solution of the following initial value 
problem of (jl.ll) 

J ut = DAu + /(it), xeM 7V ,t>-n, 
| u(x,-n) = p n (z), x G R N . 

Then, by Lemma 12.51 we have 

u(x, t) < U n {x, t) < K for all x eR N andt> -n. 

The following result provides the appropriate upper estimate of U n (x,t). 

Lemma 2.12 Assume (Ao)-(A3). The function U n (x,t) satisfies 

U n (x,t) < min{K,n(x,t)} for all x G R N and t > -n, 

where H(x,t) is defined in Theorem \2.1(A 

Proof. Let v + (x,t) = min {K, n(x, i) }. From Proposition 12.21 and Lemma l2.11l we have 

v + (x,—n) = min |K, n(a;, — n)| 

I 

= min {K, Xiv(\i(ci))e x ^ x ^- c * n+h '1 + X l+iv* : e **H»+*«+i)} 
i=l 

> ip n (x) = U n (x,-n), \/xeR N . 



By Lemma I2.5f ii). it is sufficient to show that v + (x,t) is a supersolution of on 
[—n, +oo), that is, 

v + (x,t)>T(t + n)v + (x,-n)+ f T(t - s)Q(v + (x, s))ds, Vx e R N ,t > -n. (2.10) 

J —ri 

Note that Q(u) = f(u)+Lu is non-decreasing in u for < u < K. For any x G W N ,t > —n, 
we have 

ft 

,+ t 



T i {t + n)vJ{x,-n)+ j T^t - s)Qi(v + (x, s))ds 

J —n 

<T l {t + n)K i + I T t (t - s)Qi{K)ds 

J —n 
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< e~ L{t+n) K t + f e~ L{t - s) LKids = K { . 

J —n 

Consequently, 

T(t + n)v + (x, -n) + I T(t - s)Q(v + (x, s))ds < K, Vx G R N , t > -n. (2.11) 

J —n 

Note also that A(0)v* = X*v* and 

A(Ai(ci))u(Ai(c()) = M(Ai(q)MAi(q)) = c i Ai(q)^(Ai(c i )), i = 1, • • • ,1. 
It is easy to see that the function U(x,t) satisfies the linear equation: 

n t = z>An + / / (o)n(a;,t). 

Then, for any x G R , t > —n, II(a;,t) satisfies the integral equation: 

n(x, t) = T(t + n)U(x, -n) + f T(t- s) [f'(0)U(x, s) + LU(x, s)) ds. 

J — n 

By the assumption (A3), we obtain 

Q(v + (x,t)) = f(v+(x,t)) + Lv + (x,t) 

< f'(0)U(x, t) + Lv + (x, t) < /'(O)n(z, t) + LU(x, t), 

and hence 

T(t + n)v + (x,-n) + I T(t - s)Q(v + (x, s))ds 
J —n 

< T(t + n)n(x,-n) + f T(t - s)[f'{0)U(x,s) + LU(x,s)]ds 

J —n 

= IL(x,t). (2.12) 

Combining (|2.1ip and ()2. 12[) . ()2.10|) holds and the assertion follows from Lemma [2.51 This 
completes the proof. □ 

Remark 2.13 We note that if f(u) < f'(0)u for u G [0,K], then Lemma [2. 121 is a direct 
consequence of Lemma I27T1 In fact, by f(u) < f'(0)u for u G [0, K], we have 

C7 t n < DAU n + f'(0)U n , Mx G R N ,t > -n. 

Noting that U n {x, —n) = ip n (x) < IL(x, —n) for all x G l w and 

lit = £>aii + /'(o)n(x, t), Vx el w ,t> -n. 

It follows from Lemma [2771 that U n (x,t) < U(x,t) and hence U n (x,t) < min {K, LT(x, t)} 
for all x G and t > —n. 
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Now we give the proofs of Theorem 12.91 and 12.101 
Proof of Theorem 12.91 By Lemmas 12.51 and 12.121 we have 

u(x,t) < U n (x,t) < U n+1 {x,t) < min{K,n(x,t)} 

for all x G and t > —n. Using the priori estimate of Lemma 12.61 and the diagonal 
extraction process, there exists a subsequence {U Uk (x, t)}^^ of {U n (x, t)} n gN such that 
U nk (x,t) converges to a function U p (x,t) = (U\ ]P (x, t), ■ ■ ■ ,U m;p (x,t)) in the sense of 
topology T. Since U n (x,t) < U n+1 (x,t) for any t > — n, we have 



lim U n (x, t) = UJx, t) for any (x, t) G W N+1 

n— >+oo 



The limit function is unique, whence all of the functions U n (x, t) converge to the function 
U p (x,t) in the sense of topology T as n — > +oo. Clearly, U p (x,t) is an entire solution of 
ffTTT]) satisfying fl2U) . 

The assertions for parts (ii)-(iii) and (vi)-(viii) are direct consequences of (|2.4p . There- 
fore, we only prove the results of parts (i), (iv) and (v). 
(i) Clearly, U p [x,t) > for all (x,t) G R N+1 . Since 

U n (x, t) > u(x, t) > u(x, -n) = <p n (x) = U n (x, -n) 

for all (x,t) G M N x [— n, +oo), by Lemma [231 we have ^U n (x,t) > for (x,t) G 
R-^ x (-n,+oo). This yields §- t U p {x,t) > for all (x,t) G R N+1 . Noting that 

= <kA(U i;p ) t + dxf^Up) (U 1;p ) t + ■■■ + d m fi(U p ) (U m . p ) t 

> diA(U i]P )t + m (U i]P )t, i = 1, ■ ■ ■ , m, 
where tuq = min <9j/j(n), we obtain for any r G M, 

i=l, ■■■ ,m,u£W 

(Ui. p ) t (x,t) > e m o(*-) f V i (x-y,t-T)(U i;p ) t (y,T)dy>Q,\JxeR N ,t>T. (2.13) 



Assume, by contradiction, that there exist io G {1, • • • , m} and (xo,io) £ M Ar+1 such that 
(Ui - p )t(xo,to) = 0, it then follows from (I2.13P that (Ui 0]P )t(xo,T) = for all r < to- 
Hence Ui 0]P (xo,t) = Ui - p (xo, to) for all t < to, which implies that iimi_ ! ._ 00 Ui - p (xo, t) = 
Ui - p (xo,to). But following from (|2,4p . 



Ui - P (x ,t ) > and lim U io - p (x ,t) = 0. 

£ — ^ — oo 

This contradiction yields that j^U p (x,t) » for all (x,t) G R N+l . 
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Next, we show that U p (x,t) < K for all (x,t) G R N+1 . Let V(x,t) = K - U p (x,t), 
then < V(x,t) < K and V t (x,t) < for all G l^ 1 and 

V t {x, t) = DAV(x, t) - f(K - V(x, t)). (2.14) 

We claim that V(x,t) 3> for all (x,t) G H. +1 . If this is not true, then there exist 
io G {1, ■ ■ ■ ,m} and (xo,to) G R^" 1 " 1 such that Vi (xo,to) = 0, and hence AVi (xo,to) > 0. 
It follows from (j!TT4"]) that 

< d io AV io (xo,t ) 

< f io (K-V(x ,t )) 

= fi {Ki - Vi(x ,t ), •• • ,K io -i - V io -i(xo,t ),Ki ,K io+ i - Vi 0+ i(x ,t ), ■ ■ ■ ,K m ) 

< A,(K) = 0, 

which is a contradiction. Thus V"(x, 3> and hence U p (x,t) <C K for all (x,i) G M Ar+1 . 
(iv) When xi+i = 1 5 by (|2.4p . we have 

max | max^Xj^Ci(x ■ + Cit + hi),T(t + j 
I 

< U p (x,t) < Y,xM^i(ci))e Xlic ' ){x - u ' +c ' t+hl) +v*e x * {t+h <+ l) . 
i=i 



Noting that 



lim r(t)e~ A *' = v* and lim $ Cl (0e~ Al(Cl)? = v(Ai(cj)), t = 1, • • • , I, 

it suffices to show that cAi(c) > A* for any c > c*. In fact, since A(A) > A(0) for any 
A > 0, M(A) > M(0) = A* (see, e.g., [29, Corrollary 4.3.2]). In view of M(Ai(c)) = cA x (c) 
and Ai(c) > for any c > c*, we obtain cAi(c) > A* for any c > c* and the assertion 
follows. The proof of part (v) is similar to that of part (iv) and omitted. This completes 
the proof of Theorem 12.91 

Proof of Theorem 12.101 (i) We only prove the case that U Po (t) converges to U Pl (t) in 
the sense of topology T as hi — > — oo, and uniformly on (x,t) G T\ a . The proofs for the 
other cases are similar. 

For (xi,--- >Xl+l) = (!,•■■ > 1 )> we denote ip n (x) by p£,(x) and U n (x,t) by U^(x,t), 
respectively. Similarly, when (xi>-* - iX/+i) = (0)1)"" >1)j we denote ip n (x) by <p pi (x) 
and U n (x,t) by U pi {x^t)^ respectively. Let 

Wfot) = U£ Q (x,t) - U£(x,t), (x,t) G R N x (-n,+oo), 

then < W^"(a;,t) < K for all (x,t) G l w x (-n,+oo). In view of /'(«) < /'(0) for all 
u G [0, K], we get 

^ = DAr + /([/; (M))-/(^(M)) 
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= DAW n + f'(u; () (x,t) + (1 - 8 3 )W n (x,t))W n (x,t) 
< DAW n + f'(0)W n (x, t), Vx G R N ,t > -n, 

where #3 G (0, 1). Define the function 

W(x,t) = v(Ai(ci))e Al(ci)( ^ 1+Cl ' +/ll) , (x,t) G R N+1 . 

Since 

A(A 1 (c 1 ))u(Ai(c 1 )) = M(A 1 (c 1 ))w(Ai(c 1 )) = ciAifaMAi (<*)), 
direct computations show that 

^ = DAW + f'(0)W(x, t), Vx G R w , t G R. 
Moreover, by Proposition 12.21 we have 

W n {x,-n) = U^ (x,-n)-U^(x,-n) 

< Q Cl (x ■ v\ — c\n + hi) 

< «(Ai(ci))e Al(ci)(:c - 1/1 - cin+/ll) = W{x, -n). 

It then follows from Lemma 12.71 that 

< W n (x,t) = U™(x,t) - U^{x,t) < W{x,t) = v {Xi{ci))e Xl ^ x - Ul+Clt+h ^ 
for all (x,t) G R^ x \—n,+oo). Since lim Ug(x,t) = U Pi (x,t), i = 0,1, we get 

< U P0 (x,t) - U Pl {x,t) < v(Xi(ci))e Xlici)(x ^ +Clt+hl) for all (x,t) G R* +1 , 

which implies that U po (x,t) converges to U pi (x,t) as hi — > —00 uniformly on (x,t) eT] 
for any A, a G R. For any sequence h\ with h\ — > —00 as ^ — > +00, the functions t/ p £ (x, t), 
Po := ( c i) h\, vi, ■ ■ ■ ,ci,hi,vi,hi + i), converge to a solution of (jl.ip (up to extraction of 
some subsequence) in the sense of topology T, which turns out to be U Pl (x,t). The limit 
does not depend on the sequence h\, whence all of the functions U po (x,t) converge to 
U pi (x,t) in the sense of topology T as hi — > —00, and the assertion of this part follows. 

The proofs of parts (ii)-(iii) are similar to that of part (i), and omitted. Moreover, the 
proof of part (iv) is straightforward. This completes the proof of Theorem 12.101 

3 Entire solutions for non-cooperative systems 

In this section, we consider the entire solutions of (jl.ip with monostable and non- 
cooperative nonlinearity. We introduce two auxiliary cooperative reaction-diffusion sys- 
tems and establish some comparison arguments for the three systems. Then, we prove the 
existence and qualitative properties of entire solutions using the comparison theorem. 

Throughout this section, in addition to (Ao) and (Ai), we also make the following 
assumptions: 
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(A 2 )' There exist K ± = (iff, • • • , K±) > with < K < K < K+ and two continuous 
and twice piecewise continuous differentiable functions / + ,/~ : [0,K + ] — >■ M m such 
that / G C 2 ([0,K+],R m ), /±(0) = /+(K+) = /"(K") = 0, and 

/"(«) < /(«) < / + W for all u G [0,K+]. 

(A3)' There is no other positive equilibrium of J 1 * 1 between and K ± , and f(u) and f^(u) 
have the same Jacobian matrix f'(0) at u = 0. 

(A 4 )' djf±(u) > for all u G [0,K+], 1 < j ^ i < m. 

(A 5 )' For any A; G Z+, pi, ■ • • ,p fe > and Ai, • • • , X k G [0,A*], 

/+( min {K+, p^(Ai) + • • • + p^(A fc )}) < /'(0) [pi«(Ai) + • • • + p fc «(A fc )] • 

Remark 3.1 Clearly, if f + (u) < f'(0)u for u G [0,K + ], then (As)' holds. We remark 
that when (jl.ip is cooperative, then = f and K ± = K. We also note that if / is 
defined on [0,+oo) m , then (As)' can be replaced by (A5)*: 

(A 5 )* For any k G Z+ pi, ■ ■ ■ , p k > and Ai,-- - , X k G [0,A*], 

/+ (pi^(Ai) + • • • + Pk v{\ k )) < /'(0) [pi^(Ai) + • • • + Pk v(X k )] . 

Denote W + = [0,K + ]. It is easy to verify that for any 93 G [0,K + ]x, system (jl.ip 
admits an unique solution u(x,t;ip) satisfying u(-,r;ip) = </?(•) and < u(x, t;tp) < K + 
for all x G 1^ and t > r. 

Now, we consider the following two auxiliary cooperative reaction-diffusion systems 

u t = DAu + f + (u), x€R N ,t£R, (3.1) 
u t = DAu + f~(u), x€R N ,t€R. (3.2) 

Take L = max u6 p^+ j =1 ... m \dif^(u)\ and define 

Q{u) = (Qx(n), • • • , 4») = /(«) + Lu, « G 

Q±H = (Qf(u), • • • , Q±(u)) = + Lu, u G . 

Clearly, Q ± (w) is non-decreasing in u for it G and 

Q~{u) < Q(u) < Q + (u) for any u G W + . 

We further define the operator T(t) = (T\(t), • • • , T m (t)) as (|2.2j) by replace L with L. 

The following comparison theorem plays an important role in the proof of our main 
result for the non-cooperative system. 
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Lemma 3.2 Let n,n ± E C(R N x [t,+oo),W + ) be such that 

u~(x,t) <f(t-T)u~(x,T) + J f(t-s)Q~(u~(x,s))ds, VxeR N ,t>r, (3.3) 

u(x, t) = T(t — t)u(x, t) + / T(t - s)Q{u(x,s))ds, VzeR JV ,t>r, (3.4) 

u + {x,t) > f{t - t)u + {x,t) + J f(t-s)Q + {u + (x,s))ds, VxeR N ,t>r, (3.5) 

and u~(x,t) < u(x,r) < u + (x,t). Then, there holds 

u~(x,i) < u(x,t) < u + (x, t) for all x E R N and t > r. 

Proof. We first prove u(x,t) < u + (x,t) for all x E and t > r. Let w(x,t) = u(x,t) — 
u + (x,t) and define 

dQf(u) 

Li = max — ^ , i = 1, • • • , m, and [r]+ = maxjr, 0} for any r£l, 

«ew+,i=i,-,m auj 

Since w(-,t) < and Q + (u) is non-decreasing in it for u E by (|3.4j) and (13.5H . we 
obtain 



Wi(x,t) < T i (t-r)w i (x,r) + J T { (t - s)[Qi(u(x, s)) - Qf(u + (x, s))]ds 
< j %{t - s) [Qt(u(x, a)) - Qf(u + (x, s))) ds 

fi(t-s)( f ^r<3+(u + (x, s) + Ow(x, s))de] ds 







= J Ti(t- s)(^2wj(x,s) -^-Qf{u + {x,s) + dw{x,s))dd^ds 
< J fi(t- s) ^L i jr[w j (x,s)]^ ds, Vx &R N ,t> t. 

Consequently, 

[wi(x,t)]+<J fi{t-s) ^L i ^[w j (x,s)] + j ds, Vx&R N ,t > r. (3.6) 

Let w{x,t) = Ys£=li w i(. x >t)]+- It follows from f)3. 6j) that 

m „t ____ 

w(x,t) < / Ti(t — s)Lim{x, s)ds 

i=i ^ r 
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/•t p 

< l}Z Li^i(x - y,t - s)w(y,s)ds 
= P(x-y,t- s)uj(y,s)ds, 

Jt JR n 



where P(y,s) = Y^Li Li^i(y, s). Using the same argument as in [30, Lemma 3.2], we 
obtain w(x,t) = 0, and hence u(x,t) < u + (x,t) for all x € R^ and t > r. Similarly, we 
can prove that u~(x,t) < u(x,t) for all x £ M. and t > r. This completes the proof. □ 

The following result is a direct consequence of Lemma 13.21 see also File [9] . 

Corollary 3.3 Let u^u^ G C(R N x [t, +oo), W + ) be such that Ui,uf is C 1 in t and C 2 
in x. If 

Ut < DAu~ + f~(u~), Vx £R N ,t> r, 
u t = DAu + f(u), Vx£R N ,t>T, 
uf > DAu + + / + (n+), Vx € R N ,t > r, 

and u~(x, r) < n(x, r) < u + (x, t) ; i/ien, 

u~(x,t) < u(x,t) < u + (x,t) for all x G R^,* > r. 

From the argument of Wang |36[ Theorem 2.1], we have the following result. 

Proposition 3.4 Let (Ao)-(Ai) and (A2)'-(A5) / hold. For any c > c* and v 6 R N with 
||f || = 1, ( 13. /ias a non- decreasing traveling wave solution 



& c (x-v + ct) = (</> ljC (x -V + Ct)," - , mjC (x • ^ + ct)) , 
which satisfies &c(') ^ 0> ^c(~ °°) = 0, 3>|T(+oo) = K~ and 

lim $"(e)e" Al(c)5 = v(Xi(c)), < w(Ai(c))e Al{c)f /or aZZ (el. (3.7) 

5— >— 00 

i/ere, c* ; Ai(c) and v(Xi(c)) are given as in Section 1. 

We also consider the following ordinary differential system 

u'(t) = f'(u), teR. (3.8) 
By Lemma 12.111 the following result holds. 

Lemma 3.5 Let (Ao)-(Ai) and (A2)'-(A5)' hold. There exists a solution T~(t) : R — > 
W + of (|3.8p which satisfies r~(— 00) = and r~(+oo) = K . Furthermore, 

-f-r~(i) > 0, lim T-(t)e~ x '" t = v* and T~(t) < e XH v* for all t G R, 
dt t->~ 00 

where A* = M(0) and u* = u(0). 
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The following theorem contains the main results of this section. 

Theorem 3.6 Let (A )-(Ai) and (A 2 )'-(A 5 )' hold. For any I G Z + , , v\ G R N with 

\\ui\\ = l,hi,--- ,hi+l GR, ci,--- ,q > c*, andxir" £ {0,1} withxi-\ hXH-i > 

1, f/iere exists an entire solution U(x,t) := [Ui(x, t), ■ ■ ■ ,U m (x,t)) of (jl.ip such that 

u~(x, t) < U{x, t) < min {K + , U(x, t)} (3.9) 

for all (x,t) G R N+ \ where 

u~~(x,t) = max \ max Xi$Z ( x ■ v% + erf + hi) , xi+iT~(t + h t+1 ) \, 

li=l,--,l 1 ' J 

n(a;,t) = J2xAXl(ci))e x ^ x - Vi+Cit+h ^ + X i+iv* e x ^ t+h '+^ . 
i=i 

Furthermore, the following statements hold: 

(i) U(x,t) > for (x,t) G R N+l and lim^-oo supiun^ ||£/(x,i)|| = for any A G R+. 

(ii) If Xl+i = 1; i/ien lim inft^. +00 inf xg R U(x, t) > K~ and for every x G R N , 

U(x,t) ~ !)*e A * (t+ '"+ l) as t -> -oo. 

(hi) If xi+i = 0, iaen liminf^ +00 inf i|a;||<i4 U(x, t) > AT - /or any A G R+ and for every 
x G 

I7(z,t) = 0(e^'-' Cl ^) as t -> -oo, 
where i)(a, - ■ ■ , q) = min (ciAi(ci), • • • , qAi(q)}. 

Proof. Let VF n (x,t) = (W±(x,t), ■ ■ ■ ,W^(x,t)) be the unique solution of the following 
initial value problem 

j ut = DAu + /(«), xGM 7V ,t>-n, 
I u(x, — n) = <$2 n (x), x G R^, 

where 

^ n (x) := max j^max^Xi^^ ■ v% - c^n + hi) , Xm r ~(-n + fy+i) j • 

We hrst show the following claim. 
Claim. The function iy n (x,i) satishes 

u~(x,t) < W n (x,t) < u + {x,t) : = min{K + ,n(x,t)} for all x G R N ,t > -n. (3.10) 

In fact, from Proposition 13.41 and Lemma 13.51 we see that 

u~(x, -n) = tp n (x) = W n (x, -n) < min{K + ,n(x, — n)} = u + (x, -n), Vx G R. 
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By Lemma [321 it suffices to show that for any x € Mr, t > —n, 



u 



u 



~(x,t) <T(t + n)u~{x,-n) + f T(t - s)Q'{u~ (x, s))ds, (3.11) 

J —n 

+ (x,t) > f (t + n)u + (x, -n) + I f(t-s)Q + (u + (x,s))ds. (3.12) 

J — n 

Now we prove ()3.11|) . Note that the function u(x, t) := Xj^c 3 { x ' v j JrC 3^ Jr ^ l j) (j ' = 1> ' ' ' > 0> 
satisfies the equation 

5* = DAu + f~(u), 

or the integral equation 

5(x, i) = T(t + n)u(x, — n) + / T(t — s)Q~ (u(x, s))ds. 

J —n 

Since u~(x,t) > u(x,t) for x € R^,t > —n, and Q~(u) = f~(u) + Lu is non-decreasing 
in u for u € W + , we have 

T(t + n)u~(x, -n) + / f (t - s)Q~(u~(x, s))ds 

J — n 

T(t + n)u(x, —n) + / T(t — s)Q~~ (u(x, s))ds 



> 



= u(x,t), Mx G M , t > -n, 
that is, 

T(t + n)u~ (x, —n) + / T(£ — s)(5~(n _ (a;, s))ds > Xj^7 { x ' v i + + fy?) • (3.13) 

J —n 

Similarly we can show that for x € M. N ,t > —n, 

f(t + n)u"(x,-n)+ f T(t-s)Q-(u~(x,s))ds>xi+iT~(t + hi +1 ). (3.14) 

J — n 

Hence, (l3TTTj) follows from (IBTTSj) and (IXTID . 

Next, we prove (|3.12p . Since Q + (u) = f + (u) + Lu is non-decreasing in u for u € 
we get for x € R-^ , i > — n, 

Ti(t + n)uf(x,-n)+ f %{t - s)Q+{u + (x,s))ds 

J — n 

< e - l{t+n ^K+ + I e- Z ^- s) K+Lds = K+, i = 1, • • • , m. 

J —n 
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Consequently, 

f(t + n)u + (x, -n) + f f(t- s)Q + (u + (x, s))ds < K+, Vx G R N , t > -n. (3.15) 

J — n 

Note that H(x, t) satisfies the integral equation: 

n(x, t) = f(t + n)U(x, -n) + f T(t - s) [f'(0)U(x, s) + LU(x, a)] ds. (3.16) 



By the assumption (As)', we obtain 

Q + (u + (x, t)) = f + {u + (x, t)) + Lu + (x, t) < /'(O)n(x, t) + LU(x, t). 
It follows from (|3.16p that 

f(t + n)u + (x,-n) + f f(t - s)Q + (u + (x,s))ds 



<T(t + n)U(x,-n) + / T(t - s)[f'{0)U(x, s) + LU(x, s)]ds 
= U(x,t). 



(3.17) 



Combining (I3.15P and (I3.17p . fj3. 12f> holds. Therefore, the claim follows from Lemma 13.21 
Moreover, W n (x, t) satisfies the regular estimates as in Lemma [2.5| that is, there exists 
a positive constant M, independent of n, such that for any x G 1^ and t > —n + 1, 



dW n 



Of 



x,t) 



d 2 W n 



dtx; 



and 



d 2 W n 



(xX, 



(x,t) 



d 3 W r 



d 2 W n 



dt 2 



dxft 



(xX, 



(x,t) 



d 3 W r 



dxi 



(x,t) 



d 2 W n 



dxd 



dx 2 Xj 



(x,t) 



<M, Vt,j = l, 



(x,t) 



< M, 



By using the diagonal extraction process, there exists a subsequence {W nk (xX;}keN of 
{W n (x, i)}neN such that W Uk (xX; converges to a function 

U{x,t) = (t/i(x,t),--- ,U m (x,t)) 

in the sense of topology T. Clearly, U(x,t) is an entire solution of (jl.ip . By virtue of 
(|3.10p . we have 

u~(x,t) < U(x,t) < min{K+,n(x,t)} for all (x,t) G R N+l . 

From (13.91) . it is easy to see that the assertion of part (i) holds. Note that cAi(c) > A* 
for any c > c*, and 



lim T~(t)e 

t— y— oo 



- XH = v\ lim *-(Oe-W c *K = v(\ 1 (*)), t = l,.-. X 



The assertions for parts (ii) and (iii) are direct consequences of (|3.9p . The proof is com- 
plete. □ 
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4 Applications 

In this section, we apply our main results developed in Sections 2 and 3 to the models 

4.1 A buffered system 

Consider the buffered system (jl.2p . For simplicity, we consider the case n = 1, i.e. 

) <9 t t>i = d 2 Avi + ki(b — vi) - k 2 uivi, 

where dx,d,2,ki,k2,b are positive constants. Our choice of the function g is the typical 
monostable nonlinearity, i.e. g(ui) = u±(l — u\). Let w\ = u\ and w 2 = b — v\, then (14. 1|) 
can be transformed to 

{d t wi = diAwi + wx(l — wi) + kiw 2 - k2W\{b - w 2 ), ^ ^ 

d t w 2 = c?2Aw 2 - kiw 2 + k 2 wi(b - iu 2 ). 

System (|4.2p has only two equilibria = (0, 0) and K = (l, A;2&/(^2 + &i)) and is cooper- 
ative on [0,K]. Let D = diag(di, cfe), and 

f(wi,u>2) = (wi(l - wi) + kiw 2 - k 2 wi(b - w 2 ), -k\W2 + k 2 wi{b - W2)) ■ 



Theorem 4.1 If di > d 2 , 1 > &2& and k% > k2, then the conclusions of Theorem \2.9\ are 
valid for ((4~2j) . 



It is easily seen that 
Obviously, /'(0) is cooperative and irreducible, and 



s (j/( )) = 1 ~ fc 2^ - fci + V(l - fc 2 b - fcx) 2 + 4fci > Q 



Hence, the conditions (Ao), (Ai)(a) and (A2) hold for (|4.2p . Moreover, for any A > 0, 

A: 

/c 2 & c^A 2 - fci 



^(A):=DA 2 + /'(0)=( dlA2 + 1 -^ 6 fcl 



Direct computation shows that 
M(A) = s{A(X)) 
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rfiA 2 + d 2 A 2 + 1 - k 2 b - h + yjg - d 2 )A 2 + 1 - k 2 b + fcx] 2 + 4fc 2 fcife 

2 

and the eigenvector v(X) corresponding to M(X) is 

v(X) := (v 1 (X),v 2 (X)) = (M(X)-d 2 X 2 + h,k 2 b) » (0,0). 



> 0, 



Take c* = inf 



M(A) 

a>o 



. Next, we check the condition (A3)* (see Remark 12. ip . Note that 



d\ > d 2 , 1 > k 2 b and for any A > 0, 
m(A) M (A) - d 2 A 2 + £4 



«2(A) 



1 



> 



2k 2 b 
1 



(di - d 2 )A 2 + 1 - k 2 b + h + V[(di - d 2 )A 2 + 1 - /e 2 6 + fci] 2 + 4A; 2 fei6 

> h. 



i-k 2 b + h + ^[l-k^ + ktf + ihhb 



2 L 

For any /c € Z + , pi, • • • , pk > and Ai, • ■ ■ , A& G [0, A*], denote 

(zx,z 2 ) := (piVi(Xi) H h /3jfc«i(Afc),/oiu 2 (Ai) H h pkV 2 (X k )) > (0,0). 

Consequently, (A3)* is equivalent to the following two inequalities 

zi(l - zi) + fciz 2 - /c 2 zi(6 - z 2 ) < (1 — k 2 b)z\ + k\z 2l 
-k\z 2 + /c 2 zi(6 - z 2 ) < k 2 bz\ - k\z 2 



or 



z\ > k 2 z 2 and — k 2 z\z 2 < 0. 



(4.3) 



Since > &i for any A > 0, we have z\jz 2 > k\. Therefore, (|4.3p holds if k\ > k 2 . 

4.2 An epidemic model 

Consider the epidemic model (11.31) . Scaling time and absorbing the appropriate con- 
stants into u 2 , system fjl .3|) can be rewritten as 



&tui(x, t) = d\Aui(x, t) — ui(x, t) + ju 2 (x, t), 
d t u 2 (x,t) = d 2 Au 2 (x,t) — (3u 2 (x,t) + g(ui(x, t)), 



(4.4) 



where d\ = di/an > 0, d 2 = d 2 ja\ x > 0, 7 = ai 2 /a 2 1 > and /3 = a 22 /ai\ > 0. For 
convenience, we denote dj by dj, i = 1, 2. 
We assume 



(Hi) 5 € C 2 ([0,+oo),[0,+oo)), 5 (0) = g(k) - &k = 0, > for u € (0,fc), and 

S'('u) < g'{0)u for it G [0, k], where k > is a constant. 
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(H2) One of the following holds: 

(a) g(u) is increasing for u > 0; 

(b) There exists a number ii max > such that g(u) is increasing for < u < u r 

and decreasing for u > u m3X . 

Let K = (k, g(k)/ f3), D = diag(di, cfo), and 

/(«!, U 2 ) = (-«! + 7U 2 , -/8U2 + • 
Clearly, /(0) = /(K) = and 



/'(0) 



-1 7 
ff'(0) -/3 



From (Hi), we see g'(0) > & > 0. It is easy to see that f(u) < f'(0)u for u G [0, K], f(0) 
is cooperative and irreducible, and 

s(/ / (0)) = + + V(/? + l) 2 + 4( 7 g'(0)-~fl) > Q 

Thus, the conditions (Ao) and (Ai)(a) hold for (j4.4p . Furthermore, for any A > 0, 

' diA 2 - 1 7 \ 



A(A) := -DA 2 + /'(0) 
and 



9 '(0) d 2 A 2 - /3 



M(A) = S (A(A)) = + * A ' ~ ' ~ 1 + ^ - S - ^' ~ »™ > 0. 

Clearly, inf / s i>0 Mip. exists and denote by c*. 

Theorem 4.2 Assume (Hi). The following statements hold: 

(i) J/(H 2 )(a) or (H 2 )(b) ZioWs and fc < u max , then the conclusions of Theorem \2.9\ are valid 
for (|4.4p . //, m addition, g'(u) < g'(0) for u G [0, then the conclusions of Theorem 
WM hold true for ggp . 

(ii) 7f (H 2 )(b) /ioWs and fc > n max; t/ien i/ie conclusions of Theorem \3.6\ hold for 



If (Hi) and (H2)(a) or (H2)(b) hold and k < u max , then system (|4.4p is cooperative on 
[0, K]. It is easy to verify that (A2)-(Aa) hold. If, in addition, g'(u) < g'(0) for u G [0,fc], 
then f'{u) < /'(0) for u G [0,K]. Therefore, the statement (i) of Theorem 14.31 holds true. 

When (Hi), (H2)(b) hold and k > n max , system (|4,4p is non-cooperative on [0, K]. 
Take 

n min = inf |n G (0, u max ] g(u) = g(^g(u max fj j. 
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Clearly, u m - m > 0. We define two functions f ± {u) as follows: 

f±(u) = ( - ui +7U2, -/3u 2 +g ± {ui)), 



where 



and 



9 + {ni) 



9 (ui) 



g(ux), ui £ [0, 

"max j 

<7("max)) "1 G ['Umax) ^(^max)] 

ui g [0, 

^min 7 

S^^min)) "l G [^min? jg S^^max)] • 



Clearly, g+( Ul ) < </(0)m for ui G [0, ^(n max )] . Hence, /+(«) < f'(0)u for u G [0,K+] 
which yields that (As)' holds. One can further check the conditions (A2)'-(A4)' with 
K=(k,g(k)/P), 

K + = (^g(u m3iX ),g(u mSiX fj and K = (J^g(u T ^),g(u mio )j. 

Therefore, the statement (ii) of Theorem 14.31 holds true. 
We remark that two specific functions 

uu uu 
gi(u) = — and g 2 {u) 



1 + UU 1 + vu z 

which have been widely used in the mathematical biology literature, satisfies the above 
conditions for a wide range of parameters u and v. In fact, we have the following state- 
ments: 

(a) if W7 > (3, then the function 



f(u 1 ,u 2 ) = ( - -ui + 7U2, -/3u 2 + gi(ui)) 

1=1 

j3u 



satisfies the conditions (Hi) and (H 2 )(a) with k — ^ 7 13 ■ 



(b) if W7 > f3, then the function 

f(ui, u 2 ) = ( - ui + 7U2, -/3u 2 + 52(^1)) 
satisfies the conditions (Hi) and (H 2 )(b) with 



fc = V^/ and Umax = 

Furthermore, it is easy to see that if W7 < 2/3, then k < u max , and if u^f > 2/3, then 

k > "max 1 
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4.3 A population model 

Consider the model (|1.4p by taking the non- monotone Ricker function u\e~ Ul as h{u\). 
Let w\ = u\ and w 2 = «2 — 1, then (11 .41) reduces to 



!d t wi = diAwi + wi(ri - a - 5wi + riw 2 ), ,^ ^, 

d t w 2 = d 2 Aw2 + r 2 (l + w 2 )[-W2 + h(w{)], 

where h{w\) = wie~ Wl and <ii, <i 2 , ri, r 2 , a, <5 are all positive parameters. Similar to [36J, 
we assume 

r\ > a, d\ > <i 2 and 5 > . (4.6) 

ri + r 2 - a 

In the nonnegative quadrant, (|4.5p has only two equilibrium = (0,0) and K = (K\,K 2 ) 
which satisfy 

ri K ie - Kl =6K 1 + a-n and K 2 = K x eT K K (4.7) 
Let D = diag(di,d 2 ) and 

f(w) = {w\(ri - a- 5wi + riw 2 ), r 2 (l + u> 2 )[-u> 2 + w\e~ Wl \). 

For any A > 0, 



A(X) :=DX 2 + f'(0) 



d\\ 2 + ri — a 

r 2 d 2 X 2 - r 2 



Direct computation shows that M(A) = d\\ 2 + r\ — a > and the eigenvector v(X) 
corresponding to M(A) is 

v{\) := (fi(A), ua(A)) = ((di - d 2 )X 2 + n + r 2 - a, r 2 ) > (0, 0). 

Hence, the conditions (Ao) and (Ai)(b) hold for (|4.5p . Take c* = inf^ >0 . Note 
that h(wi) = w\e~ Wl achieves its maximum at h m = 1, and is increasing on [0, h m ] and 
decreasing on [/i m ,+oo). 

Theorem 4.3 Assume ()4.6|) . The following statements hold: 

(i) //i^i < 1, then the conclusions of Theorem \2.9\ are valid for |^.5[ ). 

(ii) If Ki > 1, then the conclusions of Theorem \3. 6\ hold true for ( f^.5[ ). 

When K\ < 1, system (|4.5|) is a cooperative system on [0,K], i.e., (A 2 ) holds. We 
need to check the condition (A3)* (see Remark l2.ip . For any k £ Z + , p\, ■ ■ ■ , pf. > and 
Ai, • • • , Afc £ [0, A*], denote 



(zi,z 2 ) := (pifi(Ai)H h PfcWi(A fc ),piD 2 (Ai) H h/9fc^2(Afc)) »(0,0). 
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Consequently, (A3)* is equivalent to the following two inequalities 



or 



Z\[r\ — ot — 8z\ + T\Z 2 \ < (n - a)zi, (4.8) 
r 2 (l + z 2 )(-z 2 +z 1 e~ Zl ) <r 2 ( Zl -z 2 ) (4.9) 



6zt > nz 2 , (4.10) 
e Zl { Zl + zl) > zi(l + z 2 ). (4.11) 



Since for any A > 0, 



vi(X) _ (di - d 2 )\ 2 + n + r 2 - a r% + r 2 - a 
v 2 {\) r 2 ~ r 2 

we have 

z\ r\ + r 2 - a 
Z2 ~ r 2 

Note also that z\ > and e Zl > 1 + z\. Thus, the following two equalities suffice to verify 
dUDD and gT[]): 

S\ +r 2 - a ^ , 2 / 1^2, 3 2 



r-2 



> n and z x z 2 + [z\ - -z 2 ) + -z 2 > 0, 



which are true provided that (|4.6p holds. 

If K\ > 1, system (14. 5p is non-cooperative on [0,K]. Similar to [20l|36], we define two 
functions f (u) as follows: 



/ ± (w) = (wi(ri - a- 5wi + riw 2 ),r 2 (l + w 2 )[-w 2 + /i ± (wi)]), 

where 

h+(wi) = 



u'le""' 1 , Wx G [0, 1], 
e _1 , > 1, 



and 

,_, v / w ie ~ w \ wi G [0,/to], 

Here if^ > K\ and /to G (0, 1] are the unique roots of the equations 

SKf + a-r x - nh + {K^) = and /i e" h ° - K^e~ K i = 0, 

respectively. It is easy to verify that (A 2 ) / -(A 4 ) / hold with K = (K\,K\e~ Kx ) and K± 
(Kf,Kfe~ K ^ ), where K{ G (0,^1) is the unique root of the equation 

SKy + a-n- n/T (iff) = 0. 
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Next, we check the condition (A5)* (see Remark 13. ip . Let 

(2:1,2:2) := (pifi(Ai)H h pkVi(Xk), pi^2(Ai) H 1- PkV 2 (\k)) >(0,0). 

Consequently, (A 5 )* is equivalent to the following two inequalities 

z\[r\ - a — 5zi + r\Z 2 ] < (n - a)zi, (4.12) 
r 2 (l + z 2 ){ - z 2 + /i + (2i)) < r 2 (zi - z 2 ). (4.13) 

Note that (|4T8|) and (|4T9|) hold and /i + (2i) = 2ie~ 21 for 21 G (0,1]. To verify the above 
two inequalities, we only need to show (|4.13p holds for z\ > 1, i.e., 

(1 + z 2 )(-z 2 + e _1 ) < 21 - 2 2 , 
e(zi + zf ) > 1 + 22 for z\ > 1. 

2(1 + 4) > 1 + -22, 



that is, 



It suffices to show that 



which holds obviously. 

5 Conclusion and discussion 

In this paper, we consider the front-like entire solutions of m-dimensional monos- 
table reaction-diffusion systems in M. N . In the cooperative case, the existence and qual- 
itative properties of entire solutions are established using comparison principle. In the 
non-cooperative case, the existence of entire solutions is proved by citing two auxiliary 
cooperative systems and establishing some comparison arguments for the three systems. 
Uniqueness and stability of entire solutions of such systems seem to be very interest- 
ing and challenging problems. Besides, the issue of entire solutions of general bistable 
reaction-diffusion systems remains an open problem. 

We mention that the assumption (d±, ■ ■ ■ ,d m ) 3> := (0, • • • , 0) G R m (i.e. (jl.ip is 
non-degenerate) is crucial for our main results. When some but not all diffusion coeffi- 
cients are zero (i.e. (jl.ip is partially degenerate), system (jl.ip has weak regularity and 
compactness. For example, if di = for some i G {1, • • • ,m}, then m is not smooth 
enough with respect to x due to zero diffusion coefficient and hence the prior estimate 
for itj is not valid (see Lemma I2.6H . Recently, in [16] , we considered the entire solution 
of the reaction-diffusion system modeling man-environment-man epidemics with bistable 
nonlinear ity: 

I ^ = -pv(x,t)+g(u(x,t)). K ' ' 

To obtain the entire solution, we established the following prior estimate of solutions of 
doTT]) . see [Ml Theorem 3.3]. 
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Proposition 5.1 Suppose that w(x,t) = (u(x,t),v(x,t)) is a solution of \5.1\) with initial 
value ip £ [0,K]x, then there exists a positive constant M > such that for any ip £ 
[0,K]x, x £ R andt > 1, 

\u t (x,t)\ < M, \u tt (x,t)\ < M, \u tx (x,t)\ < M, \u x (x,t)\ < M, 

\u xt (x,t)\ < M, \u xx (x,t)\ < M, \u xxx (x,t)\ < M, \u xxt {x,t)\ < M, 

\vt(x,t)\ < M, \v x {x,t)\ < M, \v tt (x,t)\ < M. 

As mention above, v (x,t) in general is not C 1 in x when v(0, •) € C(IR; [O;^])- Hence, 
the estimates for v x , vt x and not valid. Here, we correct this mistake. We shall 

prove that v , i>t and u xx possess a property which is similar to a global Lipschitz condition 
with respect to x. In fact, we have the following result. 

Proposition 5.2 Suppose that w(x,t) = (u(x,t),v(x,t)) is a solution of \5.1\) with initial 
value ip = ((^1,932) £ C(M, [0,K]) ; then there exists a positive constant M > 0, indepen- 
dent of ip, such that for any x £ R and t > 1, 

\u t (x,t)\, \u u (x,t)\, \u tx (x,t)\, \u x (x,t)\ < M, 

\u xt (x,t)\, \u xx (x,t)\, \u xxt (x,t)\ < M, 
\v t (x,t)\, \v tt (x,t)\ < M. 

If, in addition, there exists a constant L' > such that for any r] > 0, sup^g^ \ip2{x + rj) — 
^2(^)1 < L'rj, then for any r\ > 0, 

sup \v(x + 77, i) — v{x, t)\ < M'r), sup \v t {x + 7], t) — v t (x, t)\ < M'rj, 

a;€R,t>l xeR,t>l 

and 

sup \u xx (x + r],t) -u xx (x,t)\ < M'r], 
x eR,t>i 

where M' > is a constant which is independent of ip and 77. 

It turns out that the results in [26] hold for the bistable partially degenerate system (|5.ip . 
More recently, we have extended the results to a class of two component monostable 
cooperative partially degenerate reaction-diffusion systems. However, it seems difficult to 
establish such results for general partially degenerate reaction-diffusion systems. Thus, 
an interesting problem is to adress the entire solutions of general partially degenerate 
reaction-diffusion systems. 
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